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Abstract 

Strong-coupling analysis of two-dimensional chiral models, extended to 15th 
order, allows for the identification of a scaling region where known continuum 
results are reproduced with great accuracy, and asymptotic scaling predictions 
are fulfilled. The properties of the large-A^ second-order phase transition are 
quantitatively investigated. 
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Recent developments in the analytical and numerical study of lattice two-dimensional 
SU(iV) X SU(A^) principal cliiral models have shown the existence of a scaling region, where 
continuum predictions for dimensionless ratios of physical quantities are substantially verified 

i- 

The scaling region begins at very small values of the correlation length: ^ ^ 2, well 
within the region of convergence of the strong-coupling series. Moreover by performing a 
nonperturbative change of variables from the standard Lagrangian (inverse) coupling (3 
to 

-I 

where E is the internal energy of the system (which can be measured in a Monte Carlo 
simulation), one can find agreement in the whole scaling region between the measured mass 
scale and the prediction for the mass gap - A parameter ratio offered by the two-loop 
continuum renormalization group supplemented by a Bethe-Ansatz evaluation [^. As a 
matter of fact, this may be thought as evidence for asymptotic scaling within the strong- 
coupling regime. 

We therefore felt motivated for a significant effort in improving our analytical knowledge 
of the strong-coupling series for principal chiral models, in order to check the accuracy of 
predictions for physical (continuum) quantities solely based on strong-coupling computa- 
tions. 

As a byproduct, strong-coupling series can be analyzed to investigate the critical behavior 
of the N = oo theory, where Monte Carlo data seem to indicate the existence of a transition 
at finite /?. 

Strong-coupling calculations for matrix models are best performed by means of the char- 
acter expansion. Even the character expansion however eventually runs into almost in- 
tractable technical difficulties, related to two concurrent phenomena: 

1. the proliferation of configurations contributing to large orders of the series, whose 
number essentially grows like that of non-backtracking random walks; 

2. the appearance of topologically nontrivial configurations corresponding to group in- 
tegrations that cannot be performed by applying the orthogonality of characters and 
the decomposition of their products into sums. 

We shall devote an extended paper to a complete presentation of our 15th-order strong- 
coupling character expansions for arbitrary U(A^) groups (extension to SU(A^) is then ob- 
tained by the techniques discussed in Refs. Q and |[l|]). We only mention that we had to 
rely on a mixed approach: problem 1 above was tackled by a computer program generating 
all non-backtracking random walks, computing their correct multiplicities, and classifying 
them according to their topology; problem 2, requiring the analytical evaluation of many 
classes of nontrivial group integrals, and extraction of their connected contribution to the 
relevant physical quantities, was solved by more conventional algebraic techniques. 

We pushed our computational techniques close to their limit; in order to extend our 
results to higher orders, a more algorithmic approach would be in order, especially for 
recognition of diagram topologies and group integration. 
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In the present letter we shall only exhibit the results we obtained by taking the large- 
limit of all the quantities we computed. In the strong-coupling domain, convergence to = 
oo is usually fast, and therefore large- results are a good illustration of a phenomenology 
that repeats itself (with small corrections) for finite values of A^. 

Starting from the standard nearest-neighbour interaction 

= -2Ar/3^ReTr{t/(a;) U\x+ij)}, (2) 
we computed the large- A^ free energy 

F^^lnJdUneM-SL) (3) 

to 18th order in the expansion in powers of /3. The resulting series is 

F^(3^ + 2(3^ + 4(3^ + 19/3® + 96/3^° + 604/3^^ + 4036/3^^ 

+ 55111^- + 5^/^- + 0(0. (4) 

The two-point Green's functions, normalized to obtain a finite large- A?^ limit, are 

Gix) = j^TT{U^{x)Um. (5) 

We computed all nontrivial two-point Green's functions up to 0(/3^^). The corresponding 
information is appropriately summarized by introducing the lattice Fourier transform 

G{p) = ^ G{x) exp{ip ■ x) (6) 

X 

and evaluating the inverse propagator 

G-\p) = Ao+p'A^ + m'-p')A2,o+p'A2,2+p'm'^p')As,,+p'As,s 

+ {{p^f -p^fA.^.+p^ip-'f -f)A,,,+fA,,,+f{{p''f -p'fA,,r 

+ -pV5,3 + ..., (7) 

where 

Ao = 1 - 4/3 + 4/32 - 4/33 + 12/3^ - 28/3^ + 52/3^ - 132/3^ + 324/3® - 908/3^ + 2020/3^° 
- 6284/3^^ + 15284/3^2 - 48940/3^^ + 116596/3^^ - 393052/3^^ + 0(/3^^), 

Ai = /3(1 + /3^ + 7/3^ + 4/3^ + 33/3^ + 32/3^ + 243/3® + 324/3^ + 1819/3^° + 2520/3^^ 
+ 14859/3^^ + 23124/3^^ + 123767/3^^) + 0(/3^^), 
>l2,o = -/3^ (1 + 6/3^ + 8/3^ + 57/3^ + 116/3^ + 500/3^ + 1152/3^ + 5173/3® + 11600/3^) 
+ 0(/3i^), 

^2,2 = -2/3® (1 + 2/3 + 12/3^ + 35/3^ + 121/3^ + 408/3^ + 1424/3^ + 4244/3^) + 0(/3^^), 
>l3,i = (1 + + 26/3^ + 144/3^ + 488/3^ + 1802/3^) + 0{f3^^), 
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^3,3 = 2/?" (l + 2/? + 20/5^ + 72/?=^ + 272/5^) + 0(/5^^), 
^4,0 = + 37/^2 + 84/3^) + 0(/5i 

A,2 = (l + 32/^2 + 64/33) + 0(/5i 
^4,4 = -2/5^^(1 + 2/5) + 0(/5i'^), 
^5,1 = 7/3^5 + 0(/3i6), 

A,,, = p'^ + 0i(3'^). (8) 

A number of physically interesting quantities can be extracted from G^^{p). In particular 
the magnetic susceptibility is 

X = Y^ G{x) = ^ = l + 4p + 12p^ + 36^3 + lOOp^ + 284/3^ + 796^^ + 2276/5^ + 6444/3^ 

+ 18572/3^ + 53292/3^° + 155500/3^^ + 451516/3^^ + 1330796/3^^ 
+ 3904924/3^^ + 11617404/3^^ + 0(/3^^), (9) 

while the second-moment definition of the correlation length leads to 



= = -^ = - - 4 + 3/3 + 2/3^-4/3^ + 12/3^-40/3^ + 84/3^- 296/3^ + 550/3^ 

- 1904/3^° + 3316 /3i^ - 15288 /3^2 ^ 28016 /J^^ + 0(/3^^). (10) 



and the corresponding wavefunction renormalization is Zq = Ai^. 
Moreover, by solving the equation 

G-\p^=tfIs,Py = 0) = 0, (11) 

we can compute the wall-wall correlation length, i.e. the true mass-gap fig] Eq- QTTI ) is 
algebraic in 

= 2(cosh/is - 1) = ^ - 4 + 3/3 + 2/3^ - 4/3^ + 10/3^ - 2Sp^ + 48/3^ - 206/3^ 
P 

+ 352/3^ - 1510/3^° + 2354/3^^ + 0(/3^2). (12) 

By solving 

{Px=ifJ'd/ V2, py=ifid/ V2) = (13) 
we find the diagonal wall-wall correlation length 

M! = 4 (^cosh ^ - l) = i - 4 + 3/3 + ^/33 - 3/3^ - 1^/3^ - ^/3^ + 0(/3^^). (14) 

In order to obtain the highest-order contribution to M^, it was necessary to compute a few 
(long-distance) Green's functions to 16th and 17th order. 
Lattice chiral models have a peak in the specific heat 
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which becomes more and more pronounced with increasing 0. Fig. I shows Monte 
Carlo data for the specific heat of SU(A^) models for = 21, 30, and U(A^) models for 
N = 15, 21. (We recall that U(A^) models at finite N should experience a phase transition 
with a pattern similar to the XY model, but its location is beyond the specific heat peak). 
With increasing N, the positions of the peaks in SU(A^) and U(A^) seems to approach the 
same value of j3, consistently with the fact that SU(A^) and U(A^) models must have the same 
large- limit. This should be considered an indication of a phase transition at = oo; a 
rough extrapolation of the C data indicates a critical coupling jSc — 0.306. Extrapolating the 
values of = 1/^g and (3e at the peak of C to A^ = oo, we obtain respectively ~ 2.8 
and P^e'^ ~ 0.220. 

The above picture is confirmed by an analysis of the large- A^ 18th order strong-coupling 
series of C, based on the method of the integral approximants 0J^. We indeed obtained 
quite stable results showing the critical behavior 

C^\P-f3,r, (16) 

with (3c — 0.3054 and a = 0.21, in agreement with the extrapolation of Monte Carlo data. 
Fig. |1] shows that simulation data of C approach, for growing A^, the strong-coupling deter- 
mination. 

In spite of the existence of a phase transition at A^ = oo, Monte Carlo data show scaling 
and asymptotic scaling (in the Pe scheme) even for j3 smaller then the peak of the specific 
heat, suggesting an effective decoupling of the modes responsible for the phase transition 
from those determining the physical continuum limit; this phenomenon motivated us to 
use the strong-coupling approach to test scaling and asymptotic scaling. In Fig. |^ we plot 
the dimensionless ratio fig/Mc vs. the correlation length C,g = ^/Mq, as obtained from our 
strong-coupling series. Notice the stability of the curve for a large region of values of and 
the good agreement (within a few per mille) with the continuum large- A^ value extrapolated 
by Monte Carlo data /i./Mg = 0.991(1) 0. 

In order to test asymptotic scaling we perform the variable change indicated in Eq.(0), 
evaluating the energy from its strong-coupling series 

^ = 1-^|^ = 1-G'((1,0)). (17) 

(cf. Eq. (^). The asymptotic scaling formula for the mass gap in the (3e scheme is, in the 
large- A^ limit, 

fi ^ y|l6exp AeAPe), AeAPe) = y^S^exp (-877/3^5); (18) 

/Se can be expressed as a strong-coupling expanded function of P by means of Eq. (|17|). 
In Fig. 1^ the strong-coupling estimates of fig/ -^e, 21 and Mg/ Ae,21 are plotted vs. (3e, for a 
region of coupling corresponding to correlation lengths 1 < ~ 3. These quantities agree 
with the exact continuum prediction within 5% in the whole region. 



Since the large- /3-function in the /3s scheme is not singular (as shown by Monte Carlo 
data [m and strong-coupling analysis) and the specific heat has a divergence at (3c-, then the 
relationship 

(iE{TE) = ^^C{T)(3,{T) (19) 

between the /3-function in the standard scheme Pl^T) and in the Pe scheme Pe(Te) leads to 
an non-analytical zero of Pl(T) at Pc- Pl(T) ~ |/3 — /^d", explaining the observed behavior 
in P of the Monte Carlo data for the mass gap at large . 

This phenomenon is further confirmed by an analysis of the strong-coupling series of x 
and Mq. Assuming they are well-behaved functions of the energy, we should have 

d^T^X din Ml; 

-Ip • (2°) 

Analyzing the corresponding series by a modified integral approximant scheme which forces 
the approximant to have a singularity at /5 ~ 0.3054, we found a critical behavior consistent 
with Eq. (0) {a ~ 0.2). 
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FIGURES 



FIG. 1. Specific heat vs. (3. The solid hne represents the resummation of the strong-coupling 
series, whose estimate of the critical (3 is indicated by the vertical dashed lines. 

FIG. 2. fis/Mc vs. = ^/Mg- The dashed line represents the continuum result from Monte 
Carlo data. 

FIG. 3. Asymptotic scaling test by using strong-coupling estimates. The dotted line represents 
the exact result ([l8|) . 
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